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Abstract 

We calculate the 0{a 2 s ) gluonic operator matrix elements for the twist-2 operators, which 
contribute to the heavy flavor Wilson coefficients in unpolarized deeply inelastic scattering 
in the region Q 2 ^> m 2 , up to the linear terms in the dimensional parameter e, (D = 4 + e). 
These quantities are required for the description of parton distribution functions in the 
variable flavor number scheme (VFNS). The 0{a 2 e) terms contribute at the level of the 
0(«g) corrections through renormalization. We also comment on additional terms, which 
have to be considered in the fixed (FFNV) and variable flavor number scheme, adopting 
the MS scheme for the running coupling constant. 



1 Present address: Instituto de Fisica Corpuscular, CSIC-Universitat de Valencia, Apartado de Correros 22085, 
E-46071 Valencia, Spain. 



1 Introduction 



Both unpolarized and polarized deep-inelastic structure functions receive contributions from light 
partons and heavy quarks. In the unpolarized case, the charm quark contribution may amount 
to 25-35% in the small x region, [1]. Since the scaling violations in case of the heavy quark 
contributions differ significantly from those of the light partons in a rather wide range starting 
from lower values of Q 2 , a detailed description of the heavy quark contributions is required. In 
the FFNS the corresponding Wilson coefficients were calculated to next-to-leading order (NLO) 
in a semi-analytic approach in [2]. 2 Consistent QCD analyzes to 3-loop order require the 
description of both the light and the heavy flavor contributions at this level to allow for an 
accurate measurement of the QCD scale Aqcd in singlet analyzes [4] 3 and the measurement of 
the parton distribution functions. The calculation of the 3-loop heavy flavor Wilson coefficients 
in the whole Q 2 region is currently not within reach. However, as noticed in [6], a very precise 
description of the heavy flavor Wilson coefficients contributing to the structure function F 2 (x, Q 2 ) 
is obtained for Q 2 > 10 TOq, disregarding the power corrections oc (mg/Q 2 ) 1 *, k > 1, which covers 
the main region for deep-inelastic physics at HERA. In this case, the Wilson coefficients are 
even obtained in analytic form. The heavy flavor Wilson coefficients factorize into universal 
massive operator matrix elements (OMEs) A^fx 2 /rrig) and the light flavor Wilson coefficients 
Cj(Q 2 /fx 2 ) [7] in this limit, 

H t (Q 2 /fi 2 ,m 2 Q /fi 2 ,z) = A Jt (m 2 Q / f i 2 ,z)^c i (Q 2 /fi 2 ,z), % = 2, L . (1) 

Here, fx denotes the factorization scale and z the longitudinal momentum fraction of the parton 
in the hadron. 

In the strict sense, only massless particles can be interpreted as partons in hard scattering 
processes since the lifetime of these quantum-fluctuations off the hadronic background ruf e oc 
l/{k\ + ttiq) has to be large against the interaction time r int oc l/Q 2 in the infinite momentum 
frame, [8]. In the massive case, 7iif e is necessarily finite and there exists a larger scale Qq below 
which any partonic description fails. From this it follows, that the heavy quark effects are 
genuinely described by the process dependent Wilson coefficients. Since parton-densities are 
process independent quantities, only those pieces out of the Wilson coefficients can be used to 
define them for heavy quarks at all. Clearly this is impossible in the region close to threshold but 
requires Q 2 /ttiq — r ^> 1, with r > 10 in case of F 2 (x, Q 2 ). For Fl(x, Q 2 ) the corresponding ratio 
even turns out to be r > 800, [6,9, 10]. Heavy flavor parton distributions can thus be constructed 
only for scales fx 2 3> nig. Their use in observables is restricted to a region, in which the power 
corrections can be safely neglected. This range may strongly depend on the observable considered 
as the examples of F 2 and Fl show. 

For processes in the high pj_ region at the LHC, in which the above conditions are fulfilled, one 
may use heavy flavor parton distributions by proceeding as follows. In the region Q 2 > IOijiq the 
heavy flavor contributions to the F 2 (x, Q 2 )-world data are very well described by the asymptotic 
representation in the FFNS. For large scales one can then form a variable flavor representation 
including one heavy flavor distribution, [11]. This process can be iterated towards the next 
heavier flavor, provided the universal representation holds and all power corrections can be safely 
neglected. One has to take special care of the fact, that the matching scale in the coupling 

2 For a fast implementation of these corrections in Mellin space see [3]. 

3 For a determination of Aqcd effectively analyzing the scaling violations of the non-singlet world data to 
0{a 4 s ), {a s = a s /(47r)), cf. Ref. [5]. 
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constant, at which the transition Nf — > Nf+1 is to be performed, often differs rather significantly 
from m Q , cf. [12], 

For the procedure outlined above, besides the quarkonic heavy flavor OMEs [6,13], the gluonic 
matrix elements are required. These have been calculated to 0(a 2 ) in Ref. [11]. Here we verify 
this calculation and extend it to the terms of 0(a 2 e), which enter the 0(a%) matrix elements 
through renormalization. The corresponding contributions for the quarkonic matrix elements 
were calculated in [14]. 4 

The paper is organized as follows. In Section 2 we summarize the relations needed to describe 
heavy flavor parton densities out of parton densities of only light flavors in terms of massive 
operator matrix elements. Furthermore, we point out terms to be added to the FFNS description 
in the MS scheme if compared to [2,6], which are of numerical relevance, cf. [16]. We also comment 
on the question of the effective number of flavors considering the renormalization of the process. 
In Section 3 the massive gluonic 2-loop operator matrix elements are presented and Section 4 
contains the conclusions. 



2 Heavy Flavor Parton Densities 

In the asymptotic region Q 2 ^> ttiq one may define heavy flavor parton densities. This is done 
under the further assumption that for the other heavy flavors the masses m Qi form a hierarchy 
m 2 Q i <C rriQ 2 <C etc. Allowing for one heavy quark of mass mg and Nf light quarks one obtains 
the following light and heavy-quark parton distribution functions in Mellin space, [11], 



f k (N f + 1, /x 2 , N) + h(N f + 1, /i 2 , N) = Al s Q [N f ,^,N). [f k (N f , /i 2 , N) + f^N f , /i 2 , TV)] 

m Q 



+A£ tQ [N,,-% r ,N].V(Nf,l',N) 



m 



(2) 



f Q (N f + l^\N) + fQ(N f + l,^,N) = A^lNf,-^^)-^^ 2 ^) 



V 2 



+A% g \Nf^,N\.G{Nf^\N) . 



m 



(3) 



Here fk(fk) denote the light quark and anti-quark densities, /q(/q) the heavy quark densities, 
and G the gluon density. The flavor singlet, non-singlet and gluon densities for (Nf + 1) flavors 
are given by 



S(iV / + l,/i 2 ,iV) = 



|NS 



S N f> „.,2 ' 



+A 



N )+ N f A ™Q\ N fi^ N 



mi 

2 



Nf,^y,N 

mi 



E(N f ,fi 2 ,N) 



4 For the first few values of the Mellin moment N the pure-singlet and non-singlet quarkonic OMEs were 
calculated to 0(a^) in Refs. [15]. 
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+ 



NfA 



N f 



p_ 

m 



2 ' 



AT 



4 fc 



A(iV / + l, y u 2 ,iV) 
G(N f + l,» 2 ,N) 



f k (N f + 1, p 2 , N) + h{N f + 1, /i 2 , N) - 



N f + 1 



■G(N f ,p 2 ,N) (4) 
E(N f + l,p 2 ,N) (5) 



A s 



N f ,-^-,N\-E(N f ,^N) 



m 



+ 4 S 



N 



Q 

l_ 

f 1 9 1 

m Q 



N] -G(N f ,p 2 ,N) . 



Here, 



4 S(PS ' S) = 0'|Of S(PS ' S) |j) = 8 i3 + ^alA ( r s(PS ' S) 



(6) 



(7) 



i=i 



are the operator matrix elements of the local twist-2 non-singlet (NS), pure singlet (PS) and 
singlet (S) operators (3^ s ( ps,s ) between on-shell partonic states \j), j — q,g and 

Mj = N f Mj . (8) 

Note that in the pure-singlet case the term 5ij in (7) is absent. The normalization of the 
quarkonic and gluonic operators obtained in the light cone expansion can be chosen arbitrarily. 
It is, however, convenient to chose the relative factor such, that the non-perturbative nucleon- 
state expectation values, S(A / /, fj, 2 , N) and G(Nf,fj, 2 ,N), obey 

Z(N f ,n 2 ,N = 2) + G(N f ,n 2 ,N = 2) = 1 (9) 

due to 4-momentum conservation. As a consequence, the OMEs fulfill the relations 

<y ^ ^ = 2) + = i,(io) 

N f Al tQ (N f ,N = 2)+A s Qg (N f ,N = 2) + A s ggtQ (N f ,N = 2) = 1.(11) 

The above scenario can be easily followed up to 2-loop order. Also here diagrams contribute 
which carry two different heavy quark flavors. At this level, the heavy degree of freedom may 
be absorbed into the coupling constant and thus being decoupled temporarily. Beginning with 
3-loop order the situation becomes more involved since there are graphs in which two different 
heavy quark flavors occur in nested topologies, i.e. the corresponding diagrams depend on the 
ratio p = m 2 /m 2 yielding power corrections in p. There is no strong hierarchy between these 
two masses. The above picture, leading to heavy flavor parton distributions whenever Q 2 ^> itlq 
will not hold anymore, since one cannot decide immediately in case of the two-flavor graphs, 
whether they belong to the c- or the 6-quark distribution. Hence, the partonic description can 
only be maintained within a certain approximation by assuming p 1. 

At this point we would like to add some remarks directed to readers who are not acquainted 
with the details of the calculation of heavy flavor Wilson coefficients. In Refs. [2,6] the calculation 
to 0(a 2 ) was performed in a scheme, in which the heavy quark insertion in the external gluon 
legs 5 were absorbed into the strong coupling constant through the relation 

a s = a s (p 2 ) [l + a s (p 2 )5a s (p 2 ,m 2 Q )] (12) 

,2 n e/2 



5a s (p 2 ,m 2 Q ) = S £ 



2(3 (N f ) 2(3 0)Q fm' Q 



+ 



l + ^e 2 + ^e 3 
8 24 



(13) 



3 In the present case, these are the terms oc Tp. 
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for Nf light and one heavy flavor. Here a s denotes the bare coupling constant a s = g 2 /(16n 2 ), 
P (N f ) = (U/3)C A -(4/3)T F N f ,Po iQ = -(4/3)7>, C A = 3,T F = 1/2 for SU(3) C . The spherical 
factor S £ = exp[(e/2)(7 S — ln(47r)] is set to one in the MS scheme. The bare coupling constant 
is thus given by 



a s (/x 2 ) 



l + a s (/T) 



2 .2(3 (N f + l) 



(14) 



If the above scheme is applied, cf. [6], the renormalized OMEs do not contain terms oc Tp. 
However, to express a s in the MS scheme, only the first term in Eq. (14) has to be used, while 
the second remains as a prefactor of the 1-loop contributions. 6 Hence, as has been lined out 
in [11] later, the latter term appears e.g. in Ag g in front of Aq^ in the case Q 2 3> rriq. This is 



encountered as well in the present paper for A^'g. Additionally, one has to do the same for the 
complete heavy flavor Wilson coefficients, leading to the extra terms 



r(l) 



m; 



(15) 



in the scattering cross section in the MS scheme. Here, 



H 



(i) 



m, 



F 2 V q 2 



u^i 7 -^-. I = ig:/>. 



Q 



2 ' 



P 



^,z] = 8T F {(3 
+ 



1 

2 



^+4z(l-z) + ' T §z(2z-l) 
1 o „ m b /_ . s .m 



Q 2 



- + + 2^(3,-1)- 4^ 



m. 



:C1 - : )i "^ 2 ln 



1 + 



1-/9 



1- 



g 2 (i-z) ' 



In 



1-/9 
l + f3 



(16) 
(17) 

(18) 



denote the leading order Wilson coefficients for massive quarks with the strong coupling constant 
taken out. In the same manner the contributions oc Tp in the non-lPI-contribution in Ag g , 
Ref. [11], have to be removed, to avoid double counting if the asymptotic representation for 
the heavy flavor Wilson coefficients is referred to. Since the light flavor Wilson coefficients are 
calculated in the MS scheme, the same scheme has to be used for the massive OMEs. It should 
also be thoroughly used for renormalization, as the case for light flavors, to derive consistent 
results in QCD analyzes of deep-inelastic scattering data. 

In Refs. [2,6,11] another contribution, which belongs to the inclusive heavy flavor contribu- 
tions to the structure functions F 2: l(x, Q 2 ), was not dealt with. To 0(a 2 ) these are heavy quark 
loop insertions on the initial state gluon line for the 1st order light flavor Wilson coefficient 
C 2(L) Q 2 )- T ne corresponding contribution is 



a^ 2 )A),Qhi 



9 

mi 



(19) 



6 The running coupling constant including heavy flavors in the MOM-scheme was presented in [17] to O(a^) 
recently, showing Applcquist-Carrazone [18] decoupling of the heavy quark contributions. 
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see also [19]. We also note that virtual corrections to A^q S , resp. H^ 2 ^ 5 , and A^'g need to be 
accounted for. In the asymptotic case Q 2 3> rriq they lead to +-functions, which regularize the 
soft singularity, cf. [6,11,13,14]. Here we always considered only one heavy quark contribution. 

The above expressions are derived for the FFNS. Charge- and mass-renormalization are 
performed multiplicatively for the observables. As evident from Eq. (13), a s (/z 2 ) has to be 
calculated for Nf + 1 flavors upon passing the Nf + 1st flavor threshold. In the FFNS the 
structure functions contain separate contributions of the strictly light and heavy flavors. The 
corresponding expressions for the Wilson coefficient contain anomalous dimensions which partly 
depend on Nf. In the case of the heavy flavor contributions to 0(a 2 ), [6, 11, 13, 14], (2-6), no 
closed light fermion lines contribute, however. The evolution of the three light flavors proceeds 
with Nf = 3. Due to this, there is no arbitrariness in the choice of Nf as sometimes anticipated 
in the literature. 



3 The Gluonic Operator Matrix Elements 

The description of heavy quark parton densities, Eqs. (2-6), requires the massive operator matrix 
elements given by the partonic on-shell expectation values (p\0 K \p), p — q,g, of the operators, 
cf. [20], 



F '° 

o v 



*S [-07^1-^2 ...D^ip]- trace terms , 
= 2r- 2 SSp[F MiaJ D M2 . . . D^F^J - trace terms . 



(20) 
(21) 



Here = <9 M — ig s t a A a ^ denotes the covariant derivative, t a are the generators of SU (3) c , ip the 
quark fields, A® the gluon fields, F^ v the gluonic field strength tensors, Sp the color trace, and 
S the operator which symmetrizes the Lorentz indices. The corresponding quarkonic operator 
matrix elements were calculated in Refs. [6,13] to 0(a 2 ) and 0(a 2 s e) in [14], respectively. 

The renormalized gluonic operator matrix elements A gq Q and A gg Q to 0(a 2 ) are given by 



l gq,Q 



= a„ 



i 99,Q 



+ (42g + 2Z m (Nf + 1) - ^ 1,(1) (^/))r- 1>(1) 



^, Q + z; g ^(N f + i)-z; g ^(Nf) 



+ Z~^\N S + 1) - Z^ 2 \N f ) + Z~W(N f + l)A g l 



+ 0(a 3 s ), 



(22) 



(i) 



+Z- 1 ' {1 \N f + l)ig> + + Z~W(N f + 1) - Z-W\N f )]r£V(N f ) 



+5a s A 



(i) 



0{a]). 



(23) 



Here Aij are the operator matrix elements after mass-renormalization has been carried out. The 
Z-factors Zij(Nf) renormalize the ultraviolet singularities of the operators and Tij(Nf) remove 
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the collinear singularities, cf. [6,11,13,14]. The terms Z g )JNf + 1) are equal to 



z; q \N f + i) 



a 

+0(a 



-JO) 

e lgq 



+ 



a 2 [I _ J0) S \ , 1 ( (0)0 , 1 (0) (0) , 1 (0) (0)A 



(24) 



l + a s 
+0(a 3 ) 



--J0) 



2 99 



~(0)~(0) , 1 (O) 2 ' 
% '99 ^ 2 99 



(25) 



In Eqs. (24,25), 7^ are the 0(a l s +1 ) anomalous dimensions and have to be taken - as well as (3 



at Nf + 1 flavors. We adopt the notation 7^ = 7j?(iV/ + 1) — ^(Nf) and define for later use 



m 



2 x e/2 

Q 



cxp 



Ck (e^ 



fc=2 



(26) 



To the operator matrix element A^q necessarily only non-lPI diagrams contribute. The un- 



renormalized OME A^q is given by 7 



2 Po>Q jo) I jgj I Q P) , ^2) 
£ 2 7 OT + 2 £ 99iQ gg>Q 



+ 0(e 2 ) . 



(27) 



The constant and O(e) contributions and a^ Q read 



^9,Q 



(2) 
2 99,Q 



4 N 2 + N + 2 



3 (TV - 1)7V(7V + 1) 
Pi 



- (2C2 + S 2 + S 2 ) 



8 8iV 3 + 13iV 2 + 27iV + 16 

9 (N- 1)N(N + 1) 2 



Si 



T (2) 



27 (N -1)N(N +1) 3 [ ' 
2 iV 2 + iV + 2 



(28) 



9(A f -i ) iv(iv + i)(- 2S '- 3S ^- s ; ,+46 - 6 ^') 



2 8iV 3 + 13iV 2 + 27iV + 16 
+7: 



9 (N — 1)N(N + 1) 

4 P 2 
81(iV-l)iV(iV+l) z 



27 (N - 1)N(N + l) 3 



with 



Pi = 43iV 4 + 105iV 3 + 224iV 2 + 23(W + 86 

P 2 = 248iV 5 + 863iV 4 + 1927iV 3 + 2582iV 2 + 182(W + 496 



(29) 



(30) 
(31) 



7 In the following we drop the overall factor [1 + (— l) N ]/2 in the operator matrix elements. 
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Here Sg = Sg(N) denote the (nested) harmonic sums, [21], 



k=i 



(32) 



The renormalized operator matrix element is given by 



A 



gq,Q 



a (o) 
Po,q7ot 



2 V / 2 V /i" 

Here the anomalous dimensions ^fq 1 ^ and 7gp 1 ^ are 
N 2 + N + 2 



a (°) 
(2) _ P0,Qlgg 

l gg,Q 9 



+ 0(a2). (33) 



1 91 



1 99 



igq 



-AC j. 
8C A 

CfTf 



(JV — 1)N(N + 1) ' 

N 2 + N + l 



Si -2 



(A/" — 1)N(N + 1)(N + 2) 



-2/3 (iV / ) , 



32 iV 2 + iV + 2 32 8N 3 + 13N 2 + 27N + 16 

:<Si 



t£> = 8C F T F 



3 (N - 1)N(N + 1) 1 9 (iV — 1)JV(JV + l) 2 

8 10 



(34) 
(35) 

(36) 



1-1 1 



16 6 4 



20 1 



+ fc A T F 



3JV-1 N N 2 N 3 N + l (N + l) 2 (N + l) 3 3N + 2_ 
23 1 19 1 2 19 1 2 23 1 



2 + 



3JV-1 3N N 2 3 N+l (N + 1) 2 3IV + 2 3 



(37) 



and = (8/3)Tp. A closer look at Eqs. (33,40) reveals, that the terms oc ( 2 cancel. The 
coefficients of the un-renormalized OME A g9j Q are given by 



i(i) 

^99,1 



i(2) 



/(e), 



(38) 



rrin \ £ 



/// 



^{7S ) 7£ ) +2/ 3o , Q ( 7 £) + 2/ 3o )} 



.(i) 

+ J9^ (2) (2) 
2e °99-Q + 



4/3 2 Q /(e) 



+ 0(e 2 ). 



(39) 



The constant and 0(e) contributions g^q and o^/q are 



tt(2) 



z (2) 



16(AT 2 + AT + 1)C 2 



A 5QN + 47 _ 
- 4 — Si 



7>CA<i 3 C 2 ^i + 3( A r_i) A r(iV+l)(Ar + 2) '27(A'>1) 
2A 



27(N - 1)N 3 (N + 1) 3 (N + 2) 
4(N 2 + N + 2) 2 C 2 



+7VC. 



(A" - 1)N 2 (N + 1) 2 (N + 2) (JV - 1)N 4 (N + 1) 4 (N + 2) 



(40) 
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^ y f A 9" 9^ 9(iV- l)iV(iV + l)(iV + 2)^ 3(iV+l) z 3(iV+l) 

4P 5 C 2 328iV 4 + 256iV 3 - 247iV 2 - 175iV + 54 

+ 9(iV - 1)N 2 (N + 1) 2 (N + 2) ~ 2 81(iV- 1)N(N + l) 2 1 

+ 



81(iV - 1)N 4 (N + 1) 4 (N + 2) 



+TfCf' 



4{N 2 + N + 2) 2 C 3 



+ 



P7C2 



3(N -1)N 2 (N + 1) 2 (N + 2) (N — 1)N 3 (N + l) 3 (iV + 2) 



4(JV- l)iV 5 (iV + l) 5 (iV + 2) 



(41) 



where 

P3 
Pa 

P 5 
P 6 

Pi 

P* 



15N 8 + 60iV 7 + 572iV 6 + 1470iV 5 + 2135iV 4 + 1794iV 3 + 722iV 2 - 24iV - 72 , (42) 
15iV 10 + 75iV 9 + 112iV 8 + 14iV 7 - 61iV 6 + 107iV 5 + 170iV 4 + 36iV 3 

-36iV 2 - 32iV - 16 , (43) 

3iV 6 + 9N 5 + 22N 4 + 29N 3 + A1N 2 + 28N + 6 , (44) 
3iV 10 + 15iV 9 + 3316iV 8 + 12778iV 7 + 22951iV 6 + 23815iV 5 + 14212iV 4 + 3556iV 3 

-30iV 2 + 288iV + 216 , (45) 

N s + AN 7 + 8iV 6 + 6iV 5 - 3iV 4 - 22iV 3 - 1(W 2 - 8N - 8 , (46) 
31iV 12 + 1867V 11 + 435iV 10 + 438iV 9 - 123iV 8 - 1170iV 7 - 1527iV 6 - 654iV 5 



+88iV 4 - 136iV 2 - 96iV - 32 . 
The renormalized operator matrix element A g9t Q reads 



(47) 



/m 



^99,Q 



a s -T F \n[ - 



Q 



+a 



(2) 

99,Q 



8 



i{ 2 ^( 7 (o) + 2/3o ) +7 (^(4ln 2 (^ 



2 V fi 2 



2/V?(7£ ) +2A))+7^7, 



1 91 >19 



0(a 3 s ) . 



(48 



We agree with the results for a^ n and d nn , , 



.(2) 



given in [11], which we presented in (28,40). 



The new terms cl^q and cSggQ, (29, 41), contribute to all OMEs through renormalization. 
With respect to the mathematical structure, o^gl gg \ q and o^^q, (28,40,29,41), belong to the 
class being observed for two-loop corrections before, [22]. In the present case even only single 
harmonic sums contribute. We checked our results for the moments N = 2, . . . ,8 using the code 
MATAD, [23]. An additional check is provided by the sum rules in Eqs. (10,11), which are fulfilled 
by the renormalized OMEs presented here and in Refs. [6,11,13]. Moreover, we observe that 
these rules are obeyed on the unrenormalized level as well, even up to 0(e), [14]. 

To describe the evolution of the parton distributions, Eqs. (2-6), the OMEs (33,48) have to 
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be supplemented by the corresponding 1PR terms 



4< 2 > _> A^+nHfi T N2 + N + 2 l^/V) (AO) 
A Q 9 - A Q9 + a s^o, Q T FN{N + 1){N + 2) In , (49) 

41q - 4V^V* 2 (^I . (50) 



m. 



Eqs. (49,50) agree with the results presented in Ref. [11]. In applying these parton densities in 
other hard scattering processes this modification also affects part of the massless hard scattering 
cross sections there, as outlined above. 



4 Conclusions 

We calculated the massive gluonic operator matrix elements A gq ^Q and A gg Q, being required in 
the description of heavy flavor parton densities at scales sufficiently above threshold, to 0(a 2 s e). 
We confirm previous results given in [11] for the constant terms and obtained newly the 0(e) 
terms which enter the 3-loop corrections to Aij via renormalization. We reminded of details of 
the charge renormalization and clarified that additional terms at 0(aj.) are to be included in the 
data analysis in the FFNS and VFNS using the MS scheme. 
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